Search for: All records

Abstract. Ice shelf calving constitutes roughly half of the total mass loss from the Antarctic ice sheet. Although much attention is paid to calving of giant tabular icebergs, these events are relatively rare. Here, we investigate the role of frontal melting and stresses at the ice shelf front in driving bending and calving on the scale of ∼100 m, perpendicular to the ice edge. We focus in particular on how buoyant underwater “feet” that protrude beyond the above-water ice cliff may cause tensile stresses at the base of the ice. Indirect and anecdotal observations of such feet at the Ross Ice Shelf front suggest that the resulting bending may be widespread and can trigger calving. We consider satellite observations together with an elastic beam model and a parameterization of wave erosion to better understand the dynamics at the ice shelf front. Our results suggest that on average frontal ablation rather consistently accounts for 20±5 m yr−1 of ice loss at Ross Ice Shelf, likely mostly due to wave erosion and smaller-scale, 𝒪(100 m), foot-induced calving. This constitutes only ∼2 % of the total frontal mass loss (since near-front ice velocities are ∼1000 m yr−1). Observational evidence suggests that sporadic larger events can skew this rate (we document one foot-induced calving event of size ∼1 km). Stresses from foot-induced bending are likely not sufficient to initiate crevassing but rather act to propagate existing crevasses. In addition, our results support recent findings by Buck (2024) that additional bending moments, likely due to temperature gradients in the ice, may play a role in driving frontal deflections. The highly variable environment, irregularity of pre-existing crevasse spacing, and complex rheology of the ice continue to pose challenges in better constraining the drivers behind the observed deformations and resulting calving rates. 

Accurately modeling the deformation of temperate glacier ice, which is at its pressure-melting temperature and contains liquid water at grain boundaries, is essential for predicting ice sheet discharge to the ocean and associated sea-level rise. Central to such modeling is Glen’s flow law, in which strain rate depends on stress raised to a power ofn= 3 to 4. In sharp contrast to this nonlinearity, we found by conducting large-scale, shear-deformation experiments that temperate ice is linear-viscous (n ≈1.0) over common ranges of liquid water content and stress expected near glacier beds and in ice-stream margins. This linearity is likely caused by diffusive pressure melting and refreezing at grain boundaries and could help to stabilize modeled responses of ice sheets to shrinkage-induced stress increases. 

Abstract Fast glacier motion is facilitated by slip at the ice-bed interface. For slip over rigid beds, areas of ice-bed separation (cavities) can exert significant control on slip dynamics. Analytic models of these systems assume that cavities instantaneously adjust to changes in slip and effective pressure forcings, but recent studies indicate transient forcings violate this—and other—underlying assumptions. To assess these incongruities, we conducted novel experiments emulating hard-bedded slip with ice-bed separation under periodic effective pressure transients. We slid an ice-ring over a sinusoidal bed while varying the applied overburden stress to emulate subglacial effective pressure cycles observed in nature and continuously recorded mechanical and geometric system responses. We observed characteristic lags and nonlinearities in system responses that were sensitive to forcing periodicity and trajectory. This gave rise to hysteresis not predicted in analytic theory, which we ascribed to a combination of geometric, thermal and rheologic processes. This framework corroborates other studies of transient glacier slip and we used it to place new constraints on transient phenomena observed in the field. Despite these divergences, average system responses converged toward model predictions, suggesting that analytic theory remains applicable for modeling longer-term behaviors of transiently forced slip with ice-bed separation. 

Glacial landforms provide a valuable record from which to study the history and dynamics of past ice sheets. Eskers record paleo subglacial hydrologic and sediment transport conditions because they are composed of sediment deposited by water flowing through subglacial channels. Despite decades of study, there is still debate about their formation mechanisms and little investigation of the differences between eskers formed over soft and hard beds. To address this complexity, we analysed eskers formed over soft beds along the southern margin of the Laurentide Ice Sheet (LIS) in the Lake Superior region. This included developing a new method to calculate the basal effective pressure gradient during esker formation along the subglacial channel using grain size estimates from a 20 m tall esker exposure. The morphometry and distribution of eskers were mapped with GIS to quantify their sinuosity and lateral spacing, and to compare those to the underlying bedrock elevation and sediment thickness. Lateral spacing decreased over time as the ice margin retreated, suggesting that melt rates increased during the LIS deglaciation. Furthermore, the relation between esker distribution and sediment thickness showed that eskers formed preferentially over thinner layers of sediment, irrespective of whether erosion occurred before their formation. The sedimentology of the Cable Esker exhibits a non‐monotonic pattern in channel boundary shear stress ranging from 10 to 300 Pa, alongside a basal effective pressure gradient fluctuating between −9 to −70 Pa m⁻¹. Negative basal effective pressure gradients are consistent with esker formation in channels close to the glacier terminus, which suggests lower water pressure than normally assumed. This, combined with dynamic water level fluctuations within the esker channel, supports the theory of the formation of eskers near the ice margin. 

Basal conditions that facilitate fast ice flow are still poorly understood and their parameterization in ice‐flow models results in high uncertainties in ice‐flow and consequent sea‐level rise projections. Direct observations of basal conditions beneath modern ice streams are limited due to the inaccessibility of the bed. One approach to understanding basal conditions is through investigating the basal landscape of ice streams and glaciers, which has been shaped by ice flow over the underlying substrate. Bedform variation together with observations of ice‐flow properties can reveal glaciological and geological conditions present during bedform formation. Here we map the subglacial landscape and identify basal conditions of Rutford Ice Stream (West Antarctica) using different visualization techniques on novel high‐resolution 3D radar data. This novel approach highlights small‐scale features and details of bedforms that would otherwise be invisible in conventional radar grids. Our data reveal bedforms of <300 m in length, surrounded by bedforms of >10 km in length. We correlate variations in bedform dimensions and spacing to different glaciological and geological factors. We find no significant correlation between local (<3 × 3 km) variations in bedform dimensions and variations in ice‐flow speed and (surface or basal) topography. We present a new model of subglacial sediment discharge, which proposes that variations in bedform dimensions are primarily driven by spatial variation in sediment properties and effective pressure. This work highlights the small‐scale spatial variability of basal conditions and its implications for basal slip. This is critical for more reliable parameterization of basal friction of ice streams in numerical models. 

Abstract Abrasion acts to smooth glacial terrains and leaves behind linear scratch-like features (striations) on bedrock landscapes. Striations are often used as measures of glacier flow directions, but their morphology can also provide information about the subglacial stress conditions that produced the features. While striations are often abundant in the field, the processes that create them can be opaque and hard to examine in situ because they occur under thick layers of flowing ice. To alleviate that difficulty and provide information for interpretation of the populations of striations that are observed in the field, we conducted a set of laboratory experiments in which a ring of temperate debris-laden ice was slid atop a planar marble bed under various contact force conditions that led to the creation of hundreds of striations. During the experiment, numerous glaciological properties were continuously measured, including the resistive drag. Following the completion of the experiments, the marble beds were extracted, and the striations were measured for length and categorized by morphological type, and a subset was measured using a high-resolution white-light profilometer. These experiments showed that, similar to field observations, type 2 striations were initially the most abundant; however, we found that type 3 striations became the most abundant at large displacements. We found good correlation between the abundance of striations as a function of displacement and measured drag as a function of displacement. When taken together, these results suggest that, in natural settings, ice flow around small roughness elements in glacier beds can “reset” the basal debris field, causing striations to become more abundant in their wake. As roughness is linked to quarrying, abrasion rates may increase in areas of increased quarrying. 

Tropical glaciers have retreated over recent decades, but whether the magnitude of this retreat exceeds the bounds of Holocene fluctuations is unclear. We measured cosmogenic beryllium-10 and carbon-14 concentrations in recently exposed bedrock at the margin of four glaciers spanning the tropical Andes to reconstruct their past extents relative to today. Nuclide concentrations are near zero in almost all samples, suggesting that these locations were never exposed during the Holocene. Our data imply that many glaciers in the tropics are probably now smaller than they have been in at least 11,700 years, making the tropics the first large region where this milestone has been documented. 

Abstract Glacier-bed characteristics that are poorly known and modeled are important in projected sea-level rise from ice-sheet changes under strong warming, especially in the Thwaites Glacier drainage of West Antarctica. Ocean warming may induce ice-shelf thinning or loss, or thinning of ice in estuarine zones, reducing backstress on grounded ice. Models indicate that, in response, more-nearly-plastic beds favor faster ice loss by causing larger flow acceleration, but more-nearly-viscous beds favor localized near-coastal thinning that could speed grounding-zone retreat into interior basins where marine-ice-sheet instability or cliff instability could develop and cause very rapid ice loss. Interpretation of available data indicates that the bed is spatially mosaicked, with both viscous and plastic regions. Flow against bedrock topography removes plastic lubricating tills, exposing bedrock that is eroded on up-glacier sides of obstacles to form moats with exposed bedrock tails extending downglacier adjacent to lee-side soft-till bedforms. Flow against topography also generates high-ice-pressure zones that prevent inflow of lubricating water over distances that scale with the obstacle size. Extending existing observations to sufficiently large regions, and developing models assimilating such data at the appropriate scale, present large, important research challenges that must be met to reliably project future forced sea-level rise. 

Ice at depth in ice-stream shear margins is thought to commonly be temperate, with interstitial meltwater that softens ice. Models that include this softening extrapolate results of a single experimental study in which ice effective viscosity decreased by a factor of ∼3 over water contents of ∼0.01–0.8%. Modeling indicates this softening by water localizes strain in shear margins and through shear heating increases meltwater at the bed, enhancing basal slip. To extend data to higher water contents, we shear lab-made ice in confined compression with a large ring-shear device. Ice rings with initial mean grain sizes of 2–4 mm are kept at the pressure-melting temperature and sheared at controlled rates with peak stresses of ∼0.06–0.20 MPa, spanning most of the estimated shear-stress range in West Antarctic shear margins. Final mean grain sizes are 8–13 mm. Water content is measured by inducing a freezing front at the ice-ring edges, tracking its movement inward with thermistors, and fitting the data with solutions of the relevant Stefan problem. Results indicate two creep regimes, below and above a water content of ∼0.6%. Comparison of effective viscosity values in secondary creep with those of tertiary creep from the earlier experimental study indicate that for water contents of 0.2–0.6%, viscosity in secondary creep is about twice as sensitive to water content than for ice sheared to tertiary creep. Above water contents of 0.6%, viscosity values in secondary creep are within 25% of those of tertiary creep, suggesting a stress-limiting mechanism at water contents greater than 0.6% that is insensitive to ice fabric development in tertiary creep. At water contents of ∼0.6–1.7%, effective viscosity is independent of water content, and ice is nearly linear-viscous. Minimization of intercrystalline stress heterogeneity by grain-scale melting and refreezing at rates that approach an upper bound as grain-boundary water films thicken might account for the two regimes.